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Abstract 

We study the dynamical response of a single semifiexible polymer chain based on the 
theory developed by Hallatschek et al. for the wormlike-chain model. The linear viscoelastic 
response under oscillatory forces acting at the two chain ends is derived analytically as a 
function of the oscillation frequency tu. We shall show that the real part J' of the complex 
compliance J = J' + iJ" in the low frequency limit w — > is consistent with the static result 
of Marko and Siggia whereas the imaginary part J" exhibits the power-law dependence w"'"^/^. 
On the other hand, these compliances decrease as uj~'^^^ for the high frequency limit w — > oo. 
These are different from those of the Rouse dynamics. A scaling argument is developed to 
understand these novel results. 

1 INTRODUCTION 



The recent experimental advances in the manipulation of single molecules, such as optical tweezers 
and atomic force microscopy together with single- molecule fluorescence [H [21 El ll| , have enabled 
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us to carry out mechanical and relaxational measurement in the nano-scale with piconewton sensi- 

rn 

O tivity [5j in both equilibrium and non-equilibrium conditions. For example, static force-extension 



measurements of stretching of a single polymer chain have been carried out O [H [71 [8j . As a non- 
equilibrium dynamics, the viscoelastic properties or the elastic and dissipative properties have also 
been studied [HI [13l [TOl [HI [12] . Such experiments have revealed more detailed properties of single 
molecules that are difficult to obtain in bulk experiments due to the average taken over molecules 
and time. Therefore, these investigations lead to better understanding of the hierarchical struc- 
ture of soft matter and the relationship between the molecular morphology and the functionality 
of biological molecules [H [H] . 

One of the characteristic features of soft matter such as polymers or membranes is that they 
often have several length scales. Even in a single polymer chain if the chain is semi-flexible, there 
are at least two length scales, i.e., the persistence length and the total chain length. In the several 
experiments of single polymer chains, the semifiexibility, i.e., the stiffness, is an important factor 



1 



P, El [7j. In fact, the wormlike-chain model, which is a model of a semiflexible polymer [15], [161 E], 
explains many experimental results considerably better than the flexible polymer chain model [17] 
particularly in the situation such as the highly-stretching limit in the force-extension measurement 
and the high wave-number limit of the dynamic structure factor, and so on [21 [6l [9], [71 [HI [19]. We 
emphasize that the rigidity effect can be enhanced in the above limits even for flexible polymers 
with a weak stiffness and that discrepancy appears between experiments and the theory based on 
a purely flexible model [19] . Therefore, investigation of the nonlinear dynamics due to the stiffness 
is necessary not only for semiflexible polymers but also for flexible polymers. 

Despite the above fact as well as their fundamental interest in the field of mesoscopic physics and 
their importance to the material and biological application, semiflexible polymer chains have not 
been studied intensively especially for the dynamics because of the strong nonlinearity contained 
in the wormlike- chain model. Most of the theoretical studies of single polymers have been made in 
the limiting cases of either very flexible polymers or rigid rods [20J. So far, computer simulations 
have been carried out for a stiff chain or a semiflexible chain [22l [23l [211 l2T] . 

Static theories of a semiflexible polymer chain are summarized as follows. Marko and Siggia 
derived the static force-extension relation based on the wormlike-chain model [iij. Other statistical 
properties, such as the distribution function of the end-to-end distance, have also been investigated 
[25l [26l [27] . Improvement of the wormlike- chain models has been proposed to examine the static 
properties [71 [28] . 

On the other hand, as mentioned above, analytical approaches to non-equilibrium dynamics 
of a semiflexible single polymer chain are limited. Some of the previous works have employed an 
approximation of linearization for the inextensibility constraint [29j. This linearization neglects 
non- uniformity of the line tension along the chain and has been applied to stretched polymers 
[m [301 [311 [32]. 

Recently, Hallatschek et al. [33l[M] have formulated the force-extension theory for the wormlike- 
chain dynamics without linearization of the inextensibility condition introducing the concept of 
tension propagation. They consider a weakly bend situation and use a kind of multi-scale per- 
turbation methods. The theory has been applied to the relaxation of an elongated chain after 
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removing an external force [351 EE] . 

Finally, it is also mentioned that, as a previous theoretical method, the scaling approach to a 
semiflexible polymer chain [Ml [3S] , which was successful for flexible chains [3Sl IMl HD] • 

In the present paper, we develop the linear viscoelastic theory of a strongly pre-stretched single 
semiflexible polymer chain. We consider the situation such that an oscillatory force in addition to a 
constant force is applied to the two end of a wormlike- chain. Based on the method by Hallatschek 
et al. we derive the analytic representation of the complex compliance and the complex modulus. 
It will be shown that the frequency dependence is quite different from that of the Rouse model 
|40j . The preliminary results have been published in Ref. [H]. We apply a scaling analysis to 
understand the physical insight of the results. 

The outline of the paper is as follows: In Section [21 we present the dynamical model of the 
wormlike-chain and the tension-propagation equation is derived based on the method by Hal- 
latschek et al. [33l In Section [3l the complex compliance and the complex modulus are 
obtained analytically. In Section [H the compliance and the modulus in the Rouse dynamics are 
given for comparison. In Section [5l the scaling approach is applied to both the weak-bending 
wormlike-chain dynamics and the Rouse dynamics. Summary and discussion are given in Section 

m 

2 WORMLIKE-CHAIN MODEL AND THE RESPONSE 
TO THE OSCILLATORY FORCE 

2.1 Dynamics of the wormlike-chain model 

The effective Hamiltonian for the wormlike- chain is given by [15] 

,9 2 

J^WLC — 7T 





d'r 


/ ds 


Jo 


ds^ 



(1) 

with the constraint 

|r'(.,t)r = l, (2) 

where t denotes the time, s is the length along the chain from one end, L is the total length and 
r{s,t) represents the conformation of the chain. The positive constant n is the bending rigidity. 
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The prime indicates the derivative with respect to s. The constraint ([2]) can be incorporated into 
the Hamiltonian as 



dr ^ 



, ' (3) 

as 

where f{s,t) is the Lagrange multipher for the constraint ([2]) and is interpreted as the hne-tension. 
By assuming the over-damped motion, the stochastic equation of motion of a chain is given by 

Cdtr{s, t) = -Kv"" + {f{s, ty{s, t))' + g{s, t) + ^(s, t), (4) 

where the friction coefficient C is a 3 x 3 matrix with the components Qj = x, y, z) and gr(s, t) 
represents the external force. The random force ^(s,t) obeys the Gaussian white statistics: 

<e.(s,t)>=o, (5) 

< UsMAs'.i') >= 2A;bT0/(s - s')6{t - t') (6) 

with ks the Boltzmann coefficient and T the absolute temperature. The equation of motion (j4j) is 
the same as that employed by Liverpool [12] . 

A remark is now in order. A stiff filament with an internal friction has been studied where the 
friction is supposed to arise from the internal conformation rearrangement of the filament with a 
finite radius [l3]. It is emphasized here that we have not introduced such an additional friction in 
eq m As described below, the constraint eq [2] produces a strong nonlinear coupling between the 
longitudinal (parallel to the external force) and the transverse components of the conformation, 
which causes an energy dissipation whose magnitude is comparable with the typical elastic energy. 

2.2 Weak bending approximation and multiple scale analysis 



Now we follow the theory developed by Hallatschek, Frey and Kroy [331 El] • They consider the 
situation such that the chain is elongated by the force / applied to the ends. The smallness 
parameter is introduced as e = ksT / [nfY^"^ . The conformation vector r{s,t) is divided into 
two components. One is parallel to the elongation direction (along the x-axis) and the other is 
perpendicular to it, i.e., r{s,t) = (s — r\\,r±). The basic approximation is the weak bending 
approximation such that r^(s,t)^ = 0(e) ^ 1. In this situation we have rf, = (l/2)(r^)^ + O(e^). 
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Hallatschek et al. [331 IM] have introduced a concept of stored excess length defined by 

p{s,t)=^-{r'^r. (7) 

Since the parallel component of the end-to-end distance is given by R\\ = L — {r\\{L) — 'rii(O)), we 
obtain the relation 

<AR\\>{t) = -[ < Ap> {s,t)ds + o{e), (8) 
Jo 

where Ai?|| and Ap indicate the deviation from some reference state and < .. > means a statistical 
average. 

The Langevin equation (jl]) is split into two equations for r||(s,t) and r^{s,t) with the scalar 
friction coefficients C|| and (± respectively. The equation of the transverse motion is given by 

C±dtri. = -Kr"\ + (/(s, ty^y + + ^± , (9) 

where the external force g and the random force $, are divided into the longitudinal and transverse 
components as g{s,t) = {g\\,g±) and $,{s,t) = respectively. Taking the first derivative 

with respect to s for the both sides of eqlH the equation of the longitudinal motion is given by 

C||9ir[| = +(C||-a)(A-5,r^)' 

- Krf - fis, t) + (/(s, t)r||)" - - ^f, . (10) 

Note that the sign in front of is minus because of the relation r = [s — r^\,r±). In these 
expressions, o(e^/^) terms and o(e^) terms are neglected in eqOand eqUHl respectively. This set 
of equations is solved by a perturbation expansion together with the multiple scale analysis by 
introducing two scaled variables Sg = s and si = e^/^s. Noting that the ratio of the relaxation rate 
of ry to that of r± is 0(e~^/^), one may apply an adiabatic approximation for ry. Furthermore, the 
local equilibrium approximation is employed such that the degrees of freedom in the length scale 
Ss is relaxed for a given constraint for the larger scale se. In this way, one obtains the following 
set of equations 



diexp{-A{q,s,t) + A{q,s,i))^ (11) 
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and 

<Ap>{s,t) = -^d^^F{s,t), (12) 
where q is the wave number representing modulations of the conformation r±{s,t) and 

F{s,t)= [ dif{s,i), (13) 







A{q, s, t) = 2g2 \^qH + F(s, t) j /a • (14) 

The quantity < Ap > (s,t) is the bulk value of < Ap > (s,t). See Ref. [33] for details. We 
consider the situation such that the polymer chain is in a steady condition under a constant force 
/o applied at the ends till t = and then another time dependent force A/(s,t) is switched on at 
t = 0, i.e., /(s, t) = /o + A/(s, t) for t > 0. 

The tangential vector at the chain ends is approximated to be parallel to the direction of the 
external force. This is justified in the weak bend limit |33J- The time- integral of the force along 
the polymer chain is given by 

F(s,t) = Fo(t) + AF(s,t), (15) 

where FQ{t) = fot and 

AF{s,t) = [ diAf{s,i). (16) 







2.3 Characteristic length and time 

By comparing three terms in ([3]), one notes that there are three characteristic lengths 

tf-'-f (18) 

(10) 

where ip is the persistence length of the chain and ^ has a meaning of the "screening" length. In a 
linear response as we study in the present paper, the constant force /o should be used for /. The 



total length of the chain L is also a characteristic length. The smallness parameter of the weak 
bending limit e can be rewritten as follows 

,=L^k^(m"\ (20) 



This indicates that the magnitude of the characteristic lengths has a definite order for e -C 1 as 



^i<^i,<L. (21) 



Hereafter we ignore the shortest one ^f. 

Comparing each term in the Langevin equation (jl]), we obtain the following characteristic times 

Ti = ^ (22) 

K 

r2 = ^ (23) 
with £ a length scale. Substituting I = ^ into eq [23], we obtain 



n = ^ • (24) 



Substituting i = £„ into eq [221 we have 



Note that this is the only characteristic time which does not contain neither / nor L. 
2.4 Tension dynamics 

In this subsection, we focus on the propagation of the line tension f{s,t) or F{s,t). Here it is 
mentioned that this concept itself can also be applied to a flexible polymer chain [44j. Combining 
eqs [m and [12], the tension propagation equation is obtained as the closed form with respect to F; 

^ ^2^, .^ r, fl -exp(-A(g,s,t)) 







^ rft exp (-A(g, s, t) + A{q, s, i)) | . (26) 



This equation is rewritten in terms of the dimensionless quantities as 

1 -exp(-i(g, s,t)) 



KdiF{s,t) 



g2 + 1 



-2g2 ^ diexp (-i(g, s, t) + i(g, s, t)) i , (27) 



where 



q = ^q, (28) 

s = e'l'sr' , (29) 
i = t/T^ . (30) 

K = 7t/( is just a numerical factor with ( = C||/C±- The total length L is now rescaled as 
L = e^^'^L^^^. The scaled functions A and F are given by 

A{q,s,i) = 2f(^qH + F{s,i)j (31) 

and 

F{s, i) = ^F{s, t) = f dil^^^p^ . (32) 
i^^i Jo Jo 

We assume that A/(s, t) is sufficiently small and apply the linearization approximation to (l27Il . 
That is, we substitute (fT5|) into (l27|l and retain the terms up to the first order with respect to AF 
so that we obtain 

/"* 



KdiAF + / dtAF(s, t - t)M(t) = 0, (33) 
Jo 

where the memory function M(t) is given by 

M{i) ^^fdq |g4e-2''^(^^+^)* - ■ (34) 

The asymptotic behavior is given by M(i) ~ i^^ with /5 = 5/4fort— >0 and /3 = 5/2 for t ^ cxd. 

Equation fl33|) is to be solved under the boundary conditions specified by AF(0, t) and AF{L, t). 
In what follows, we consider the symmetric case that AF{0,t) = AF{L,t) = AF{t). Applying 
the Laplace transformation with respect to t to eq[33l we obtain 

KdjAF{s, z) + N{z)AF{s, z) = , (35) 

where -F(s, z) denotes the Laplace transform of -F(s, t) and 

N(z)=a[ dql^r—X—, -^^\ (36) 

Jo ^Wiq^ + l) + z 2g2 + l/ ^ ' 
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is the Laplace transform of M{t). The asymptotic form of N{z) is given as follows. For uj —>■ oo, 
from eq [36] we obtain the following equations after some manipulation 

ReN(±icjr5) = -2Si{ujTi:y/\ 

lmN{±iujT^) = T^S2{ujT^y/^ (37) 

with Si = j^dq{Aq^ + 1)-^ ^ 0.863 and S2 = dqq\Aq^ + 1)'^ ^ 0.179. It is readily shown 
that the Taylor expansion of N(z) with respect to z breaks down and therefore N(z) is not analytic 
at z = 0. The correct expansion is obtained after some manipulation as follows 

N{±iujT^) = -2S3{uJT^f^'^ T iSiiuj^y , (38) 

where Ss = 7i/8 and S4 = 7r/4. 

We consider the case that the force Af(t) at the boundaries is oscillatory as A/(t) = /a sin(u;t) 
with the amplitude and the frequency u. The scaled form of AF at the boundaries is given by 

AF{i) = iujT^)-'^ [1 - cos (curg)] (39) 
Jo 



and the Laplace transform is 



/O Z[Z^ + [UJT^y] 



The solution of eq [35] can be represented as 



cos (b{z){2s-L) 

AF{s, z) = AF{z) X (41) 

cos I B(z)L 



where B(z) = {N{z) / AKy/'^ . One needs to evaluate the inverse Laplace transform of eq[lT 

1 t-c+ioo _ cos [B{z) {2s — L)] 

AF{s,t) = — AF{z) ^ ::^e^'dz. (42) 

^7r« Jc-ioo cos I B{z)L 

This will be carried out in the next section. 

3 ANALYTICAL RESULTS 



Now, we study the response of the end-to-end distance to the oscillatory force. The average end- 
to-end distance AR{t) which is a deviation from that of the steady state under the constant force 
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/o is given by 

AR{t) = - [ ds<Ap>{s, t) (43) 

= ^{dsF{s,t)l=L-dsF{s,t)l=o} ■ (44) 
Cii 

Hereafter, for abbreviation, we represent the statistical average of the end-to-end distance as AR{t) 
without the brackets < ■ >, the bar ^ and the parallel mark -n. Substituting the solution (l42l) into 
eq HH together with eq [32l we obtain the time evolution of the average end-to-end distance under 
the given boundary condition. 

Since we are concerned with the asymptotic behavior t +00, we consider only the poles on 
the imaginary axis z = 0, ±iujT^ to carry out the inverse Laplace transform. The final result can 
be written as 

Ai?(t) _ /a 



L fo 

The scaled complex compliance is given by 



J'{uj) sm{ut) — J"{uj) cos(co't) . (45) 



2D ( 

J' (a;) = Im B(icjr5) tan(B(icjrg)) 1 , (46) 



ID ( 

J"(uj) = Re B{iujT^) tan(B(icjrg)) | , (47) 



where B(z) = aN^zY^"^ /2 with complementary dimensionless constants 
and 

1 knT _ J_ 

The scaled elastic modulus G' and the scaled loss modulus G" are obtained from J' and J" as 



D = A,^^ = -^e . (49) 



follows 



G'(^) = , (50) 

J'iuY + J"(a;)2 

G'\uj) = . (51) 

J'(ujY + J"(uY 



in 



In the following, we introduce another characteristic time. The linearized eq [33] reduces to the 
following simple diffusion equation by employing the Markov approximation; 

KdlAFis,i)-^diF = 0. (52) 

This implies that we may define a new relaxation time by the following form as the time scale of 
the slowest mode, just as the Rouse time in the continuous Rouse dynamics; 



a 



2 



By this relation, the three parameters r, and a are not independent of each other. In Figures 
[1] and [2], we choose a and r as the independent parameters. 

We examine the limiting behavior of J' and J". For the high frequency limit, substituting eq 
1371 into eqs HE] and HT] and after some manipulation, we obtain 



^1/2^1/2^5/8^7/8^ 



^ ^ <'\a{ksTY'y ^_,f, ^ «:-5/8(^^2^)+i/2^-7/8 (54) 

7^1/2^1/2^5/8^7/8^ 

and 

J\uj)/J"{uj) = b/a ^ 0.199 , (55) 

where a ~ 0.721 and b ~ 0.142 are the positive solutions of (a + bi)'^ = 1/2 + 182/81. It should be 
noted that the unsealed complex compliance J = LJ/fo depends on neither L nor /q. 

For the low frequency limit, substituting eq[38]into eqs |46] and HT] and after some manipulation, 
we obtain 

A-) - (56) 
J"ico) ~ M;|k!a;+V2 n%kBT)+'u^'/' . (57) 

Note that fISB]) is consistent with the result of Marko and Siggia for the static stress-strain relation 
[6], which is given by 

R(fo) 1 keT , , 

-y^ = ^-o^- (58) 
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From this, we have 

Rim + 6)) i?(/o) _6 knT ^ ^^^^ 



L L 

Equations ( H6|) and ( H71) give us the complex comphance as a function of u. Figures [U^a) and 
[U^b) show the comphances J' and J" for a = 1 and a = 100, respectively. As mentioned above, the 
compliances exhibit the fractional power law behavior for the high frequency and J' is consistent 
with the static result of the wormlike-chain for u; — > 0. The difference for the simple Maxwell-like 
elasticity is more evident for G' and G" as plotted for a = 1 in Figure [2] (a) and for a = 100 in 
Figure [2](b). Note that both G' and G" increase as cj''/^ for cut^ ^ 1. 

In addition, an intermediate region exists only if r ^ or a ^ 1. From the definition 
f HHj) . this condition is realized in the situation that the total chain length L is much larger than 
e^^'^ip = ^I'^^^J'^ . When this condition is satisfied, there is a finite interval of the intermediate region; 
l/r^tu^l/r^. For example, in both Figure [U^b) and Figure [2]^b), the interval 1 ^ tur < 10^ 
corresponds to this region. In this region, the asymptotic form of is given by eq[3Sl Moreover, 
since B{iujT^) = aN^iuT^Y^'^ /2 ~ (zco'r)^/^, the imaginary part of B is very large. Therefore, we 
can approximate tan(i?) by +i and, substituting eq[38]into eqs |16] and HH the compliance becomes 



Thus, the compliance has the a;"^/^ dependence in the intermediate region. 

4 COMPARISON WITH ROUSE DYNAMICS 



(60) 



In this section, following the paper by Khatri and McLeish [ID], we present the complex compliance 
for the Rouse model and compare it with the present result. The Rouse dynamics without internal 
friction is governed in the continuum limit by 

^dr(n,t) ,d'^r(n,t) . . 

= ^^b^ + ^^""^ + ' ^^^^ 

where C, is the friction coefficient and the argument n indicates the ra-th monomer from one end, 
r{n) is the position vector of the ra-th monomer and k is the elastic coefficient of the linear spring 
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Figure 1: J' and J" as a function of lot for = 1 and (a) a = 1.0 and (b) a = 100.0. The full 
curve represents J' whereas the broken curve represents J". The characteristic time r is defined 
by eq[ 
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Figure 2: G' and G" as a function of ujt for D = 1 and (a) a = 1.0 and (b) a = 100.0. The full 
curve represents G' whereas the broken curve represents G". 

between a pair of adjacent two monomers. It is noted that the argument n and the number of 
monomer are treated as real numbers and satisfy ^ n ^ N. Over-damped and Markov motion 
is assumed. Both end points are subjected to the external forces which have the same amplitude 
but the opposite direction 

f{n,t) = f{t)[S{n-N)-5{n)] . (62) 

The last term ^„ in eq [61] is the White Gaussian noise that satisfies the fluctuation dissipation 
relation of the second kind 



< ^{n, t)|T(m, t') >= 2kBTCIS{n - m)5{t - t') 



(63) 
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where / is the unit matrix and two adjacent matrices mean a tensor product. 

We define the end-to-end distance as R{t) = r{N,t) — r{0,t) and the deviation as AR(t) = 
R{t) — R{0). In the same way, the deviation of the external force is defined by A/. The complex 
compliance Jr{uj) = Jr{^^) + iJ^i^) is defined through the relation 

< AR > (icu) = Jl{uj)Af{iuj) , (64) 

where the asterisk * means the complex conjugate and Jjii^j) is given by [IQ] 

where is the Rouse relaxation time defined by 

rR = ^. (66) 

The function (!65l) is plotted in Figure [3] and the corresponding complex modulus Gr is plotted in 
Figure HI 

From the expression fl65l) . the asymptotic behavior is derived to compare with that of the 
weak-bending wormlike- chain dynamics. For u; — > oo, the complex compliance behaves as 

J';,{uj) OC , (68) 

and for 



J'fiiuj) const. (69) 
4(cj) OC cj+i . (70) 

These exponents are distinctly different from these obtained in the previous section, —7/8 in both 
and as ^ oo and +1/2 in J" as — > 0. See eqs 15^ and 1571 Moreover, the viscoelastic 
behavior of the Rouse dynamics with internal friction is also examined by Khatri and McLeish 
|4Uj . where the high frequency behavior is given by oc uj~'^ and Jj^ oc uj~^. These are again 
different from the present results. 
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10-^10-^0° 10^ 10^ 10^ 

Figure 3: The complex compliance Jr{uj) = Jr{uj) + iJ'/^^u) for the Rouse dynamics without 
internal friction. The full curve represents Jjj whereas the broken curve represents J^. The 
amplitude is scaled such that = 1 for uj 0. 
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10-2 10-^10° 10^ 10^ 10^ 

Figure 4: The complex modulus Gr{oj) = G'j^{u)+iG'^{u) for the Rouse dynamics without internal 
friction. The full curve represents whereas the broken curve represents G'^. The amplitude is 
scaled such that G"^ = 1 for cij ^ 0. 

5 SCALING APPROACH 

5.1 Scaling form of AR{t) 

In this section, we apply the scaling analysis in order to explain the behavior of complex compli- 
ances for both high and low frequency limits. 

All the parameters are scaled out in eqs [27] and [221 Therefore, the parameters appear only 
through eq [32] and through the boundary condition which contains L. The scaled form of L is 
given by L = e^^^LC,"^. These facts together with eq [44] give us the following scaling property of 
AR; 

AR{t) = e'"r\fC'Q{^"'Li-\ t/r^) , (71) 
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where Q{x,y) is an unknown function whose asymptotic form is to be determined. This scahng 
form (!7T!) is very crucial to investigate the asymptotic behavior of the comphance and the modulus 
as shown below. 

5.2 Complex compliance for the high frequency limit 

The exponent 7/8 exhibited by J' and J" for c<j ^ oo obtained in eq EH can be understood by 
the following scaling analysis. In the linear response regime, the dimensionless function Q in eq 
[71] should be proportional to Ja/ fo- At the high frequency limit, the effect of the external force is 
expected to be localized near the two ends and hence the compliance J should not depend on L. 
Therefore, from eq[7l|, the compliance takes the following form 

/^JM ~ lA,y^c\fC\u;n)-' (72) 

with an unknown exponent z. Substituting the definitions of C,, e and given, respectively, by 
([19]), (EOj) and (El) into eq [72] yields 

We can require that the compliance is independent of the screening length ^ (and hence /) in the 
high frequency limit. This is because the relaxation of the chain has the factor nq'^ + /g^ as can be 
seen from eq[l]or eqs [H and [TS] and hence k is relevant for the high frequency (/g^ is irrelevant). 
This requirement gives us 

The scaling analysis in the Rouse dynamics is different from the above because of the absence 
of the local length scale, i.e., the persistence length ip. The Rouse model has only one length scale, 
the root of the mean square end-to-end distance a. When no external force is present, it is given 
by 

/u T\ 

^ iVV2 . (75) 



(7 _ 

k 

Therefore the dimensional analysis tells us that the deviation of the end-to-end distance should 
obey 

AR ~ a^J{ujTn)^^' (76) 
Kb! 
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with the external force 

fit) = fAc''^' . (77) 

Assuming that J has a power law behavior J{ujt^ ~ (cjt/j)^^ as c<j ^ oo. The complex compliance 
becomes 

/^J(cu) ~ a^iujTR)-^ ~ fAN''^^uj-^k^-\-^ . (78) 

In the high frequency limit, the response is localized and the compliance should be independent 
of N so that the exponent is determined uniquely as z = 1/2 or 

J{uj) oc a;~^/2 . (79) 

This is the argument given by Khatri et al. [40j. 

5.3 Complex compliance for the low frequency limit 

The exponent 1/2 exhibited by J" for c<j ^ obtained in eq[57]can be understood as follows. In 
the low frequency limit, the effect of the external force is extended almost uniformly to the whole 
chain. Therefore, we can require that J is proportional to L so that 

IaAuj) = ^-fe"\-h'"Lr\fC-\ur,Y . (80) 

The real part J' is independent of the frequency for ^ and hence 2; = whereas the imaginary 
part J" should be independent of n for ^ 0. As mentioned above, the relaxation of the chain 
has the factor nq^ + fq^ and hence n is irrelevant for the low frequency. Therefore, substituting 
the definitions of S,, e and given, respectively, by (fT9l) . ( j20l) and (!24l) into eqlHOl it is found that 
the exponent is given by 

= \ («) 

so that J" oc kBTuj^/"^. Note that, from eq[73l J" at the high frequency limit is proportional to 
(kBTy^"^ whereas it is proportional to ksT at the low frequency limit. 

In contrast, the complex compliance in the Rouse dynamics is analytic in the — limit. This 
fact is clear from the expression fl65|) . This should be compared with that of the wormlike- chain 
dynamics (l56l) and (1571) where the function N(z) in B{z) contains a non-analyticity as eq[38l 
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6 SUMMARY AND DISCUSSION 



In summary, we have developed the analytical theory of the viscoelasticity of single semiflexible 
polymer chains and have obtained the linear compliance which has a frequency dependence char- 
acteristic to the semiflexible chain. In particular, it is found that the asymptotic behavior of the 
compliance obeys as J', J" oc uj^"^/^ for u; — oo whereas J" oc uj^^^"^ for ujt <^ 1. These are 
distinctly different from the results of the Rouse dynamics. The constant of for u = given by 
eq El] is also different from that of the flexible chain. 

The theory assumes weakness of the bending parameter e <^ 1 which guarantees the scale 
separation. This is due to the fact that the characteristic length parallel to the stretched chain 
i\\ ~ As and the characteristic wave length ~ i± satisfy 

^ ~ ^ ~ ^'^'^ « 1 > (82) 
£|l As As ^ ^ 

where the scaling forms (|28l) and (!29l) have been used. It is emphasized that, for u oo, the scale 

separation is valid without assuming the smallness of e because of the fact that £[[ oc tu^^^^ ':$> i± oc 

^-1/4 rjj_^jg ^^^^ jg verified by eqs 15^ and 

Now we discuss the relation between the present results and those obtained by Hallatschek et 
al. who have considered the relaxation of the end-to-end distance after step-wise change of the 
external force [Ml [36] . They have predicted that both in the stretching case and in the release 
case the end-to-end distance behaves as 

< Ai?||(t) >oc //€-^/^(A;bT)-i/V/^ , (83) 

where t ^ t/ ~ (nf^^. The exponent 7/8 is the same as that in the high-frequency limit in eq 
EH At a short interval after the force change, its effect is localized near the chain ends which is 
small compared with both the persistence length f|T7|) and the screening length f|T9|) . In fact, we 
can show that eq [83] is consistent with our eq [511 as follows. In the linear response theory the 
relaxation function ilj{t) and response function 0(t) are related to each other as (j){t) = —dilj{t)/dt. 
The complex compliance is the Fourier-Laplace transform of the response function and hence 
J{uj) = "ipiO) + iuj dte^^^ipif)- Therefore, we have the following relation J{ijj) ~ iIj{1/uj) and 
J(cj) = fJ/L ~ (Ai?||(l/cj))/L. 
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In the intermediate time region ti ^ t ^ tj with the crossover time ti defined through 
^\\(tL) = L, Hallatschek et ah have obtained [31] 



< Ai?,| it) >oc fl^K-"\kBTf'^t"^ (84) 

for a pulhng situation and 

< Ai?i| (t) >oc f'^K-^'\kBTf/h^/^ (85) 

for a release situation. We have no results corresponding to eq [H] since this contains a nonlinear 
effect of the applied force. On the other hand, the exponent 1/2 in eq [531 corresponds to eqEDlin 
the present paper. Actually one can verify that not only the exponent but also the coefficient in eq 
[551 is consistent with our result. This implies that the expression fISS]) is free from the nonlinearity 
between the force-strain relation. 

Finally we mention a theoretical study which gives us the exponent 1/2 in the compliance. 
Caspi et al. have investigated the mean square displacement of a single monomer of a prestressed 
semifiexible network [3^. They have obtained 

<A/.^(x,t)>cx^^tV^ (86) 

where h{x, t) denotes the undulation amplitude, u the line tension, rj the solvent viscosity and L the 
total chain length. Equation fl86l) holds in the time region Atttik/i''^ <^ t ^ rjL'^/u. Furthermore, 
they have shown that the effective time dependent friction Ce(^) satisfies the generalized Einstein 
relation 

ksT < Ah\t) > 

W)~ 2t • 

Combining eqs [86] and Wf\ one obtains the complex compliance (J oc t/(eit)) 

Jicu) oc z/-i/2^-i/2cu-i/2 . (88) 

Some experiments of semifiexible networks support the exponent 1/2 [171 US]. It is mentioned, 
however, that the physical of this result is different from our present result for semi-fiexible chain 
given by eqlSDl This fact is clear because the coefficient in eq[HH]does not contain ksT whereas 
our expression eq[60]is proportional to {ksTY^'^. 
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Now we comment on the several effects which have not been considered in the present paper. 
The hydrodynamic effect has not been investigated quantitatively in a nonlinear wormlike- chain 
although it is expected to be not so strong in a strongly stretched semi-ffexible chain. The previous 
studies of the hydrodynamic effect in the linearized wormlike- chain dynamics |18l [T9[ [31], [32] should 
be extended to apply to the present theory. The internal friction considered in the Rouse dynamics 
[10] should also be extended to the semi-flexible chains. In addition, the helical wormlike- chain 
model, which contains the torsional energy, has been studied in dilute solutions [26|, [521 [53] . This 
torsional effect may affect the viscoelastic properties of single polymer chains. 

Before closing this article, we make an estimation of the characteristic times and r defined by 
eqs [21] and [53] respectively. The data for A-DNA in an aqueous solution are as follows [IS] [H] [31] 

ip ~ 50 [nm], 
L ~ 20 [/im], 

C± ~ 1.3 X 10"^[Pa s] = 1.3 X 10"^[pN-s//im2] . 

For these values together with ^ ~ 1/2 for a rigid rod [20j and the room temperature ksT ~ 
4.1[pN • nm], and for the external force /o ~ 10 [pN] the characteristic times are given by 

~ 2.7 X 10-^[s] 
r ~ 6.0 X 10"^ [s] = 60 [/xs] 

and the constant a ~ 5.3 x 10^. We expect that the frequency of the order of 60 [/is] is accessible 
by atomic force microscopy and that the present predictions can be detected experimentally. 
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